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Motivated by analogue models of classical and quantum field theory in curved spacetimes and 
their recent experimental realizations, we consider wave scattering processes of dispersive fields 
exhibiting two extra degrees of freedom. In particular, we investigate how standard superradiant 
scattering processes are affected by subluminal or superluminal modifications of the dispersion re- 
lation. We analyze simple 1+1-dimensional toy-models based on fourth-order corrections to the 
standard second order wave equation and show that low-frequency waves impinging on generic scat- 
tering potentials can be amplified during the process. In specific cases, by assuming a simple step 
potential, we determine quantitatively the deviations in the amplification spectrum that arise due 
to dispersion, and demonstrate that the amplification can be further enhanced due to the presence 
of extra degrees of freedom. We also consider dispersive scattering processes in which the medium 
where the scattering takes place is moving with respect to the observer and show that superradiance 
can also be manifest in such situations. 
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I. INTRODUCTION AND MOTIVATION 

Experimental realizations of analogue black holes [TJ [5] 
and their associated effects have drawn a lot of atten- 
tion in the past few years [MS]- Probably the most dis- 
cussed results have been the first observation of the clas- 
sical analogue of Hawking radiation in an open channel 
flow [5] and the still controversial observation of radi- 
ation in ultrashort laser pulse filaments |U [BJ-fH] ■ Even 
though it remains an open question whether or not real 
black holes emit Hawking radiation, calculations involv- 
ing analogue black holes suggest that the emission pro- 
cess is most probably unaffected by transplanckian effects 
that could, in principle, alter or even exclude the radia- 
tion process [TOfT^] . 

Superradiance p!3Hl5] is another typical phenomenon 
of black hole physics [THl HZ] which is also manifest in 
analogue models of gravity [TB]. In standard scattering 
processes, the ratio between the reflected and the inci- 
dent particle number currents [19] (i.e. the reflection co- 
efficient) is smaller then one. This is directly encoded 
by the fact that the amplitude of the reflected wave is 
usually smaller than the amplitude of the incident one in 
non-dispersive normalized scattering processes. However, 
in some special situations (e.g. wave scattering in a Kerr 
black hole spacetime), low- frequency incident waves can 
be amplified in the scattering process. This amplification 
effect, known as superradiance, was first discovered by 
Zel'dovich in the context of electromagnetic waves [20] . 
and later shown to be a more general phenomenon in 
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physics [21 [T5], in which classical as well as quantum 
field excitations can be amplified. 

The main goal of our work is to analyze the robustness 
of the amplification process for dispersive fields. The dis- 
persive fields considered in this paper exhibit two extra 
degrees of freedom, such that one can study the overall 
robustness of the amplification process for multiple su- 
perradiant scattering. Additionally, in some simplified 
cases we also investigate the specific deviations in the 
amplification spectrum due to dispersion. The examples 
discussed in this paper are motivated by analogue mod- 
els of gravity, where fields exhibiting sub or superluminal 
dispersion relations arise naturally. 

A. Superradiant systems 

Before starting our analysis of superradiance for dis- 
persive fields, it is important to review the basic ingredi- 
ents that make the phenomenon possible in some simple 
systems. For example, in a Kerr black hole, incident 
scalar field modes of frequency uj and azimuthal number 
772 ^ are known to be superradiant for sufficiently low 
frequencies. Although the Klein-Gordon equation is very 
complicated when expressed in standard Boyer-Lindquist 
coordinates, the radial part of the equation of motion for 
these modes can be written, after a change of variables, 
very simply as 

!^ + K,, m (r> = 0, (I) 

where u is related to the radial part of the separated 
scalar field and r* is a tortoise- like coordinate (which goes 
to — oo at the event horizon and reduces to the standard 
radial coordinate r at spatial infinity). Kj >m (r*) is the ef- 
fective potential and possesses the following asymptotics, 
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V u , m (r m ) -> 



(lu - mtt h ) 2 , r* 



— > — oo, 
+00, 



(2) 



where fih is the angular velocity of the black hole at the 
event horizon. It can be shown that, for < lu < mfl^, 
incident modes are superradiantly scattered by the black 
hole [16l [21] • The corresponding reflection coefficient is 
given by 



\R\' = 1- 



— mQh 



lir > 1. 



(3) 



Two factors are responsible for the occurrence of the ef- 
fect [15]: first, the event horizon behaves as a one-way 
membrane that allows only ingoing transmitted waves 
(defined by the group velocity) near the black hole; sec- 
ond, because of the ergoregion of the rotating black hole, 
low-frequency ingoing waves are associated with an out- 
going particle number current at the horizon (the current 
being proportional to lu — mQh) |19j . This outgoing flux 
at the horizon is compensated by the superradiantly re- 
flected modes in such a way that the total particle num- 
ber current is conserved during the process. 

Another simple system exhibiting superradiance is a 
massless scalar field <fi with electric charge e and min- 
imally coupled to an electromagnetic potential — 
(V(x),0) in 1+1-dimcnsions [13 . The evolution of this 
system is determined by the Klein-Gordon equation, 



(d IM -ieA fl )(d> l -ieA> i )<f> = 0, 



(4) 



which, after separation of the temporal dependence 
(i.e. <f> = f(x) exp(— iLut)), reduces to 



d 2 / 
dx 2 



+ (lu - eV(x)) 2 f = 0, 



(5) 



where lu is the frequency of the mode in question. If the 
electromagnetic potential V(x) has the following asymp- 
totic behaviour, 



V(x) 



0, x — > —00, 
e$0; x — > +00, 



(6) 



one can calculate the relation between the reflection and 
transmission coefficients for incident waves 1131: 



\Rfl = l- "- e *° \T\ 2 . 



(7) 



Therefore, low frequency (0 < lu < e$o|^| right-moving 
modes originating from — 00 interact with the scattering 



It is important to remark that the electric potential is defined 
only up to an arbitrary constant C. Therefore, if we redefine 
V{x) — > V(x) + C, the supcrradiant condition will be given by 
< ui — C < e<I>o. Throughout this paper, we fix the potential 
by requiring that V(— 00) = C = 0. Consequently, the condition 
for superradiance, which is based on the smallness of u> — C, is 
reduced to a condition on the smallness of u. 



potential, resulting in transmitted right-moving modes 
and superradiantly reflected (|i?| 2 > 1) left-moving 
modes. The important ingredients for superradiance are 
essentially the same as in the rotating black hole case: 
firstly, there is no left-moving mode at +00 (this is sim- 
ilar to the boundary condition at the event horizon of a 
black hole, where no outgoing solutions are allowed); sec- 
ondly, the form of the electromagnetic potential allows for 
low-frequency right-moving modes to be associated with 
left-moving particle number currents at +00. 



II. MODIFIED DISPERSION RELATIONS 

In order to study how a modified dispersion relation af- 
fects superradiance, we shall first generalize the simple, 
non-dispersive, 1+1-dimensional model above |13) by in- 
cluding fourth-order terms in equation Q. Inspired by 
the quartic dispersion relation fl 2 = fc 2 ±fc 4 /A 2 , where A 
is a dispersive momentum scale and the ± notates super 
and subluminal dispersion respectively, we propose the 
following generalization of equation Q , 



1 *4 
=F O 



+ (d fl -ieA^)(d^-ieA^)4> = 0, 



(8) 



After separation of the temporal dependence, instead of 
equation ((Sb, one has 



=F— f" 



+ J" + (lu- eV{x)f f = 0, 



(9) 



where / represents the field mode with frequency lu > 
and prime denotes derivative with respect to x. We 
choose the effective frequency fi(x) = lu— eV(x) to satisfy 
asymptotic relations similar to those in equations ([2| and 

n 2 ( \ J for x -> -00, 

I (to — e$o) 7 tor x ~ * +00, 

where e^o is a positive constant. In the asymptotic re- 
gions, the solutions of equation ^ are simple exponen- 
tials, exp (ikx), whose wavenumbers k satisfy the disper- 
sion relations below, 



0. 



T fc 4 /A 2 - k 2 + (lu - e$ ) 2 = 0, 



x 

X 



-00, (11) 

+00. (12) 



In the subluminal case (lower sign), real solutions 
(lu, k) to the dispersion relation correspond to the in- 
tersections of a lemniscate (figure- eight) and a straight 
line, see FIG. [I] In the superluminal case (upper sign), 
real solutions correspond to the intersections of a quartic 
curve and a straight line, see FIG. |4j In order to relate 
the asymptotic solutions at +00 and —00 without hav- 
ing to solve the differential equation for all values of x, 
one needs a conserved quantity analogous to the Wron- 
skian for second order wave equations. Since our model 
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is non-dissipative, we expect such a quantity to exist [55] . 
After some trivial algebra, it is possible to show that the 
expression 



where 



Z[f] =W 1 + W 2 T^ 2 W 3 , 



Wx = f"*f-rf", 

W 2 = f*f"-f"*f, 
W 3 = f*f-f*f, 



(13) 



(14) 



generalizes the notion of the Wronskiarj^Jto our dispersive 
model, i.e. dZ/dx = for any solution / of equation 
Furthermore, it is convenient to work with the scaled 
functional X whose action on a function / is defined by 



X[f] 



Z[f] 
2iA 2 ' 



(15) 



In particular, the action of X on a linear combination of 
'on shell' plane waves (wavenumbers satisfying the dis- 
persion relation) is simply 



X 



Ak n x 



(16) 



We would like to emphasize the simplicity of the alge- 
braic expression above, which only depends on the am- 
plitudes, effective frequencies and group velocities of the 
various degrees of freedom participating in the scatter- 
ing process. This result becomes even more surprising 
when one considers the intricate form of the ad-hoc in- 
troduced generalized Wronskian ( 13 ) which is used to 
derive equation (16). Note that equation (16) remains 



valid in the non-dispersive limit (A — » oo) and, there- 
fore, we rightfully consider this quantity to represent the 
generalization of the particle number current. 




FIG. 1. The dispersion relation for the subluminal case shown 
in terms of the dimensionless variables k/A and fi/A, where 
k is the wavenumber and Q — lj — eV(x) is the effective 
frequency. The green dash-dotted line represents £1 — cu 
when x — > — oo while the blue dashed line corresponds to 
£7 = ui — e^o when x — > oo. Here < u) < e"I>o is the fixed 
lab frequency satisfying the condition for superradiance. We 
have indicated the incident mode by a black square dot, the 
transmitted and reflected modes by black circular dots and 
the remaining solutions of the dispersion relation by white 
circular dots. 



A. Subluminal scattering 

Based on the general notion of superradiance described 
in the introduction, we will study the scattering process 
of incident waves originating from x — > — oo in the pres- 
ence of a subluminal dispersion. In realistic scenarios, we 
do not expect the dispersion relation fi 2 = k 2 — k 4 /A 2 
to have a maximum/minimum value above/below which 
only imaginary solutions of the dispersion relation are 
possible. Therefore, in order to guarantee that at least 
one real mode is present in the dispersion relation, we 
assume that the dispersion parameter A is large com- 
pared to the electromagnetic interaction, i.e. we assume 



2 For the standard 2nd order wave equation, Wi = W2 = and 
the functional W3 is the conserved quantity commonly referred 
to as Wronskian. For Klein-Gordon fields, the Wronskian can 
be interpreted as the particle number current. 



that e4>o < A/2. In such situations, the behaviour of the 
solutions to equation (|9| in the asymptotic limits is cap- 
tured in FIG. [I] where the green dash-dotted (f2 = ui) 
and the blue dashed (CI = oj — e$o) lines represent the 
effective frequency when x — > — 00 and x — ¥ +00, respec- 
tively. Note that there exist four roots, corresponding to 
four propagating degrees of freedom as one would expect 
from a fourth order differential equation. We also note 
that, when uj < e$0) the blue dashed line is located be- 
low the f2 = axis and the roots in the region x — > +00 
inherit a relative sign between their group and phase ve- 
locities with respect to the roots at x — > +00 (compare 
the intersections of the green dash-dotted line (x — > —00) 
and the blue dashed line (x —> +00) with the red solid 
curve) . 

We label the four roots associated with a frequency 
< cj < e<I>o as k a < kb < k c < kd when x — > —00 and as 
< ks < kc < ko when x — > +00 (they correspond to 
the intersections of the straight lines with the figure-eight 
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in FIG. [I]). The following table summarizes the character 
of these modes, where u and v notate right-movers and 
left-movers respectively, 



Roots at x — > 


a 


b 


c 


d 




Roots at x — > 


A 


B 


C 


D 


— oo 












+00 










group velocity 


u 


V 


u 


V 




group velocity 


V 


u 


V 


u 


phase velocity 


V 


V 


u 


u 




phase velocity 


V 


V 


u 


u 


effective 


+ 


+ 


+ 


+ 




effective 










frequency 












frequency 











Since, in this setup, there is only one source of waves, 
located at x — > —00, and no incoming signal from x — > 
+00, we impose the boundary condition that the modes 
A and C are unpopulated in the scattering process. Fur- 
thermore, we choose the incoming mode at x — > —00 to 
be entirely composed of low- momentum c modes with no 
high-momentum a mode component]^] 

The corresponding solution of equation ([9]) in the 
asymptotic limits is 



e'^™" + R 1 e llc '-i -f. jx 2 



R ik r 



X — > —00, 

X — > +OO, 



(17) 



where the wavenumber of the incident mode is given by 
kin = h, the reflected modes are given by k ri = kb, 
k r2 — kd, and the transmitted modes are k tl — fee, k t2 — 



kn- In addition, the coefficients Ri, i?2, T\ and Ti can 
be related to reflection and transmission coefficients (see 
equation (19) below). 

In the asymptotic regions, it is possible to solve exactly 
the dispersion relation and find the following explicit ex- 
pressions for the roots, 



A / /4a; 2 

TiVW 1 - A 2 - 



(18) 



1-1/1 



4(cj - e$ ) 2 
A 2 ' 



ku = 



In order to compare the particle number current of 
the reflected waves with the incident and transmitted 



currents, we substitute equation (17) into the expression 
for the conserved generalized Wronskian, equation (15), 
and find, after straightforward algebraic manipulation, 
the following relation between the coefficients R\ , R2 , T\ 
and T2, 



l^i I 



|i? 2 | 2 



A 2 - 4(w - e$ ) 
A 2 - Auj 2 



fcti 



\Ti\ 



11 1 n 



|T 2 | 2 



> 1. 



(19) 



This relation should be compared to the standard re- 
sult for non-dispersive ID scattering, |i?| 2 = 1 — |T| 2 , 
and its generalization in the presence of an external po- 
tential, |i?| 2 = 1 - (ui - e$ )|T| 2 /w, see equation Q. As 
discussed above, in the non-dispcrsive case it is possible 
to achieve \R\ > 1 for sufficiently low frequency scattering 
with < lu < e$o- For fields with subluminal dispersion 
relations, the conclusion is similar. From expression ( 19 ) 



above, which is valid only when < ui < e$o < A/2, we 
conclude that the total reflection coefficient (i.e. the ratio 
between the total reflected particle number current and 
the incident current) is given by the LHS of equation ( 19 ) 
and is always greater than one, characterizing a gener- 
alization of the usual superradiant scattering which in- 
volves extra degrees of freedom. Note that, in the general 
case, without an exact solution we have no information 



3 Note that here, and everywhere else in this paper, we consider 
only stimulated scattering processes in which the ingoing high 
momentum channels are suppressed. This is not applicable to 
spontaneous scattering processes, where the quantum vacuum 
naturally supplies all ingoing, low and high, momentum modes. 



about how this total reflection is distributed between the 
low and high wavenumber channels represented, respec- 
tively, by |-Ri| 2 and \k r2 / k in \\R2\ 2 ■ However, by looking 
at the A series expansion of the reflected and transmitted 
wavenumbers, we can draw some interesting conclusions 
about the regime A 3> 1. Using equation (18), we obtain 



the following expansions for the relevant wavenumbers, 

+ 0(A- 3 ), 
O (A- 2 ) , 



h . 


— hi — 




= A -2A 




hi 


= u> — e$o 


ki 2 





CO 

2A 2 



(20) 



(lu - e$ ) J 

2A 2 
^ ^2 



2A 



C(A- 



As mentioned before, all possible modes of the sys- 
tem are characterized by the same conserved frequency 
lu. The momentum of each mode, on the other hand, is 
determined by the wavenumber k and, therefore, is not 
the same for every mode. From the expansions above, 
we note that h 2 ,k r2 ~ 0(A) while h n , k ri , h ± ~ C(A°). 



5 




co / e§ co / e§ 



FIG. 2. The reflection coefficient |-Ri| 2 for the subluminal 
case as a function of u)/e&o calculated for the step potential 
V(x) — e«l>oO(2;) with e$o = 1- Note the singular behaviour 
at cj = e$o/2, which is present even in the non-dispersive 
case. This divergence is directly related to the discontinuity 
in V(x) and would be cured by smoothing the step potential 
at x = 0. 



FIG. 3. The total reflection coefficient ji?i| 2 + \k r2 / k in \\R 2 \ 2 
as a function of aj/e<E>o for a subluminal dispersion rela- 
tion. Calculations were performed assuming a step potential 
V(x) = e$o6(i) with e<I?o = 1- Note the enhancement effect 
in the amplification caused by dispersion. 



In particular, the difference between the high momentum 



modes k T2 and kt 2 is k r 



k t2 ~ ©(A" 1 ). Intuitively, 
this suggests that the creation of a pair of these modes 
{k T2 , k t2 ) should be favoured since it requires a negligible 
momentum change in the system if A 3> 1. 

We can lend weight to the arguments above by solv- 
ing the idealized problem in which the potential is cho- 
sen to be a step function, i.e. V(x) — e$o6(a0 (see 
the appendix for a discussion concerning the appropri- 
ate boundary conditions used at the discontinuity point 
x = 0). In such a case, one can show that the reflection 
and transmission coefficients are given by, 

\Kr 2 ^r 1 }{^r 1 K t 1 )[ K r 1 K t 2 ) 
(kin ~ kri){kin ~~ ^ti)(^m ~ &t 2 ) (21) 

(k ri k T2 )(k T2 fcj 1 )(fc r2 fc* 2 ) 
i /.• . - k. Vfc,_ - /.-.. ••:/.-.. - /.-. ^ 

Ti = 
T 2 = 

where the wavenumbers k are explicitly given in equa- 
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tion (18 1. Note that the expressions above are valid for 
any value of A, not only in the regime A > 1. In or- 
der to observe the effects of the dispersive parameter, 
we plot the reflection coefficient |i?i| 2 as function of lu 
for different values of A, see FIG. [2j It is interesting 
to note that as we increase the dispersive effects (i.e. we 
lower A), the reflection coefficient associated with k ri de- 
creases. We can also see in FIG. [2] that the R\ channel 
even becomes non-superradiant for some frequencies in 
the range < ui < e^o- Of course, if we also include 
the other reflection channel, see FIG. [3j then the total 
reflection coefficient is always superradiant, as proven in 
equation (19). Note that the presence of extra scatter- 



ing channels even enhances the amplification process in 
comparison with the non-dispersive case. This can be un- 
derstood in the sense that the amplifier is more effective 
if the number of 'accessible' channels increases. 

Note also that the total reflection coefficient in the dis- 
persive case is not continuously connected with the non- 
dispersive regime. More precisely, as A is increased in 
FIG.[3j the total reflect ion coefficient becomes larger and 
larger and moves away from the non-dispersive coefficient 
(A = co). A possible explanation for this behaviour re- 
sides in the fact that the subluminal dispersion relation 
itself is not continuously connected to the linear disper- 



(> 




wave originating from x — > — oo, which is converted into 
a reflected left-moving mode, a transmitted right-moving 
mode and exponentially decaying modes (as a boundary 
condition, we impose the fact that there can be no ex- 
ponentially growing modes in the asymptotic regions). 
Therefore, the solution of equation Q corresponding to 
this scattering process is given by 



k/A 



FIG. 4. The dispersion relation for the superluminal case 
shown in terms of the dimensionless variables k/A and f2/A. 
We have chosen here a frequency in the superradiant inter- 
val < ui < e4>o- The blue dashed line corresponds to the 
dispersion relation at x — > +00 while the green dash-dotted 
line corresponds to the dispersion at x — > — 00. The incident 
mode is indicated by a black square dot, the transmitted and 
reflected modes by black circular dots and the remaining so- 
lutions of the dispersion relation by white circular dots. 



sion relation in the limit A — > 00. In other words, the 
modes k r2 and k t2 which are always absent in the non- 
dispersive case, will be present in the dispersive regime 
no matter how large A is. As argued before, the pro- 
duction of these modes is then favoured when A ^> 1, 
explaining the observed discontinuity. 



Re 



E r e k 



Te iktX + E t e~ k ' tX , 



x — > —00, 
x +00, 



(22) 



where the coefficients R and T are related, respectively, 
to reflection and transmission coefficients (see equa- 



tion (24) below), and the coefficients E r and E t are the 
coefficients of the exponentially decaying modes. The 
wavenumbers k in equation ( 22 ) are obtained directly 



from the dispersion relations, equations (II) and (12) 



71' 



-1 + 4/1 



4lu 2 
A 2 " 



k P r — 



A / 4w 2 

7IV 1 + V 1+ a^ 



(23) 



-1 + 4/1 



4(w - e$ ) 2 

A 2 



A^ 1 + «. 



The conservation of particle number current can be 



obtained by inserting equation (22) into equation (16) 



and equating the generalized Wronskian at ±00. Com- 
pared to the non-dispersive result, one might expect ex- 
tra terms related to the exponentially decaying modes. 
However, these extra contributions are also exponentially 
decaying and, at ±00, their particle number currents are 
fully absorbed in the low momentum channels. As a con- 
sequence, one obtains the following reflection coefficient 
(valid for < u < e$ ), 



B. Superluminal scattering 

Let us now turn our attention to the case of a superlu- 
minal dispersion. Once again, we can understand much 
of the scattering process from the dispersion relation, see 
FIG. 21 The main difference with the subluminal case is 
that now there are only two real roots in either asymp- 
totic regions, corresponding to one left-moving and one 
right-moving mode. The other two solutions to the fourth 
order equation are imaginary roots, corresponding to ex- 
ponentially decaying and exponentially growing modes. 
Similarly to the subluminal case, we are interested in the 
scattering of an incident low-frequency (0 < uj < e$o) 



\Rf 



A 2 + 4(w-e$ ) 2 
A 2 + 4w 2 



k {n 



\T\ 2 > 1, (24) 



which, similarly to the subluminal case, is always greater 
than one. Note also that the expression above reduces to 
the usual non-dispersive reflection coefficient |7]) in the 
limit A — » 00. 

Repeating the analysis used in the subluminal case, 
one can solve exactly the problem for a step potential 
given by V(x) = e$o@(^)- Using appropriate boundary 
conditions at x = (see the appendix) , one can relate the 
coefficients in equation ( 22 ) to the wavenumbers given in 
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equation (23), 





co/e<j) 



0.6 



0.8 



FIG. 5. The reflection coefficient |i?| 2 for the superluminal 
case as a function of u)/e&o calculated for the step potential 
V(x) = e"I>o0(a;) with e"I>o = 1- Like in the subluminal case, 
the divergence at oj — e$o/2 is caused by the discontinuity in 
V(x) and would be absent if a smooth potential were used. 
Note the enhancement in amplification due to dispersion. 
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(25) 



+ k t ) 



Once again, it is useful to plot the reflection coefficient 
as a function of ui for different values of the dispersive 
parameter A, see FIG. [5] 

The results obtained for the scattering of sub and 
superluminal dispersive fields in our 1+1-dimensional 
toy-model demonstrate that superradiance is possible 
and that the amplification is enhanced due to the extra 
degrees of freedom in comparison with non-dispersive 
superradiance. These results are interesting for experi- 
mental attempts of detecting superradiance, since they 
indicate that dispersion may increase the amplification 
rates and, consequently, make the effect easier to observe. 

III. GENERALIZATION TO NON-ZERO FLOWS 

Having in mind moving media in analogue models of 
gravity, we extend the ideas of the previous section by in- 
cluding a position-dependent flow velocity W(x) in our 
model. In other words, the medium where the scattering 
process takes place (fluid frame) is in relative motion with 
respect to the observer (lab frame). If we require that the 
dispersion relation be unaltered in the comoving frame of 
the fluid, the modified Klein-Gordon equation ^ gener- 
alizes to 



T-j^d^ + + (d t - ieA t + d x W) (9* - ieA* - Wd x ) = 0, (26) 

where the derivative operator d x is understood to act on everything to its right. After separation of variables, 
4>(t,x) — e~ lult f(x), we obtain the following wave equation, 



which can be written as 

f""(x) + a(x)f"(x) + /3(x)f'(x) + j(x)f = 0, (28) 
where 

a = ±A 2 (l - W 2 {x)) 

/3 = =fA 2 [2W(x)d x W(x) - 2iW(x) (lu - eV{x))] (29) 
7 = ±A 2 [(cj - eV(x)) 2 + id x (W(x)(u) - eV(x)))] . 



Observe that these coefficients are not independent but 
satisfy the following relations, 



7 - 7 * = 2£Xm( 7 ) = ilm(d x (3) = d x [ ^— ^ ) , (30) 



Ke(J3) = -(p + 0*) = 9 x a. 



(31) 



We choose the functions W(x) and V(x) to be asymp- 
totically constant, 



(32) 
(33) 




so that, at ±00, any solution to equation (271 can be 
decomposed into plane waves satisfying the dispersion 
relation below, 



k 2 + — - 

k ± A 2 - 



(u> - e$ - kWoY 



x — > —00, 
x — » +00, 



(34) 



where the ± stands for super and subluminal respec- 
tively. As previously, we work under the assumption that 
w,e$ < A. 



The generalized Wronskian ( 13 1 is not conserved in 
this case. However, one can show that a modification of 
it, namely 



Z = W 1 + W 2 +aW 3 -ilm(/3)|/| : 



(35) 



and the corresponding scaled quantity X — Zj (2iA 2 ) are 
independent of x. The action of the functional A on a 
linear combination of 'on-shell' plane waves takes pre- 
cisely the same form as in the zero-flow case, see equa- 
tion ( fl6] ) (the only difference is that the effective fre- 
quency is now given by Q = oj — eV — kW instead of 
fl = oj — eV). This result highlights the generality of 
the algebraic structure of the particle number currents 
given by ( [l6| . Additionally, the existence of superradi- 
ance in the zero-flow case was related to the condition 
that oj — e$o < 0. We can anticipate that the occurrence 
of superradiance in non-zero flows will be favoured by 
modes for which oj — e$o — kWo < 0. 



A. Subluminal dispersion 

Let us consider first the case of a subluminal disper- 
sion relation. As shown in FIG. [6] (red solid curve), for 
fixed e<I>o/A < Wo < 1, there are two distinct intervals of 
frequencies separated by a critical frequency w cr it < e$ 
in which we expect superradiance: for < oj < w cr j t (re- 
gion I) two propagating degrees of freedom are admitted, 
both right-moving in the lab frame; for w cr it < oj < ei>o 
(region II) there are four real roots of the dispersion re- 
lation corresponding to four propagating degrees of free- 
dom, three right-moving and one left-moving (with re- 
spect to the lab-frame). Note that the requirement that 
W is not too small, specifically W > e$ /A, is nec- 
essary in order to guarantee that the left-most root in 
region I (i.e. the circular dot labeled k t2 in region I of 
FIG. [6]) has positive group velocity and hence defines a 
true transmitted mode. If < W) < e<E>o/A, the situa- 
tion is basically the same as the one discussed in section 



II A On the other hand, if Wo > 1 (see the black dashed 
curve in FIG. |6|, there is only one possible regime: for 
all frequencies < oj < e<I>o, two right-moving modes are 
admitted. 
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FIG. 6. The subluminal dispersion curve in the lab frame 
(A = 1, e$o = 0.4) at x — > +00 for two different flow ve- 
locities: e$o/A < Wo < 1 (red solid curve) and Wo > 1 
(black dashed curve). The coloured regions described in the 
text are associated with the red solid curve: region I (green, 
< oj < LJcrit), region II (blue, u; cr it < oj < e«I>o), and region 
III (light red, non-superradiantly scattering region). The in- 
tersections of the horizontal lines with the red solid curve and 
the black dashed curve indicate transmitted modes at that 
frequency. Note that, when x — > — 00, the dispersion is de- 
scribed by the green dash-dotted curve of FIG. [l] 



An interesting fact is that, due to the absence of left- 
moving modes at +00 when oj < u; cr i t , this system is a 
model for the event horizon of an analogue black hole 
with modified dispersion relations. However, the pre- 
cise location of the horizon, besides being w-dependent, 
is also rather ill-defined, relying on a global solution to 
the equation of motion in the vicinity of a classical turn- 
ing point. This region can be studied by WKB methods 
and Hamilton- Jacobi theory [23] , but it is not of specific 
interest to us here. 

According to the analysis above, the scattering of 
an incoming wave from —00 will result in transmission 
through two or three channels, depending on whether 
oj < aj cr it or not. An exact solution to the scattering 
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problem can be decomposed as 



r ie ifc « x + T 2 e ik * 2X , 



R 2 e l 



x -t 



-oo, 

-00, 



(36) 



when < lo < u; cri t (if Wo < 1) or < u < e$ (if 
Wo > 1) and as 



Rte 



ik r ix 



R 2 e 



k r 2X 



X 



-oo, 



Tie 



iktix 



X 



-oo, 



(37) 



.T 2 e ikax +T 3 e iktsX , 
when uj cr it < u < e$o (only possible if Wo < 1 



The wavenumbers fa n , k ri and k r% are given by equa- 
tion ( 18 1 and are labeled in FIG.[l] while the transmitted 
wavenumbers are not, in general, expressible as simple 
functions of the parameters. Note that we do not explic- 
itly keep track of possible exponential decaying solutions 
in equation ( 36 ) (corresponding to complex solutions to 
the dispersion relation), since they do not contribute di- 
rectly to p8l). 



We insert the solutions (36) and (37 1 into the func- 



tional X of equation ( 16 1 and find, after algebraic ma- 



nipulations, the following relationship between the trans- 
mission and reflection coefficients, 



\Ri\* = l- 



A 



VA 2 - 4w 2 



\ n Km 



u-e$ -fa n W )\T n \ 



(38) 



where the sum is over all (2 or 3, depending on lo and Wo) 
transmission channels. Here, v 9n — v g (fa n ) are the group 
velocities of the transmitted modes at +00. The LHS of 



equation (38) can be interpreted as the total reflection 



coefficient associated with the incident modes. 

Since the group velocities v 9n of the transmitted modes 
are, by definition, always positive, the sign of the con- 
tribution from each channel fa n to the RHS of (38) is 
determined by the factor fi(fc) = to — e$o — kWo eval- 
uated at kt n , as we anticipated previously. When only 
two transmission channels are admitted, one of the two 
factors f2(fc) is strictly negative while in the case of three 
transmission channels, two of the three factors 0(fc) are 
strictly negative. Therefore, in both situations there is 
one root which contributes an overall negative amount 
to the RHS of equation ( 38 1 and thus reduces the mag- 
nitude of the total reflection. Because of these trouble- 
some modes, the RHS is not strictly greater than 1 and 
we cannot straightforwardly conclude superradiance. In 
general, in order to fully answer the question of superra- 
diance, one would need to specify W(x) and V(x) in all 
space and solve for the coefficients R n and T n . 

Large A approximation: To better understand the 



relation between equation (38) and superradiance, we 



now focus on small deviations (A » 1) from the non- 
dispersive case. For fixed Wo < 1, the size of re- 
gion I in FIG. [6] becomes zero when A is greater then 
e$ [2(1 - Wo)/3]~ 2/3 . We therefore assume that the fre- 
quency to lies in region II, where three transmission and 
two reflection channels are available, see equation (37). 



As explained above, in the most general case one would 
need to solve the equation of motion for all x in order to 
conclude superradiance or not from equation ( 38 ) . 



Similarly to the zero flow case, by looking at the series 
expansions of the relevant wavenumbers, we can make 
useful predictions about the transmission coefficients in 



I 

this A 3> 1 regime. We start the analysis by solving 
equation ( 34 ) in the asymptotic region x — > +00 and ex- 



pressing the obtained transmitted wavenumbers as power 
series in A, 



fa, = 



LQ-e^o 1 (^-e$ ) 3 A _ 2 
1 + Wo 2(1 + Wo) 4 



W< 



to — e$ 
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1 



wl 



-0(A- 4 ) : 
O (A- 1 ) 



(39) 
(40) 



The series expansions of the incident wavenumber fa n and 
of the reflected wavenumbers, k ri and fc r2 , are given, as 
previously, by equation ( |20[ ). From these A expansions, 
we note that, unless Wo ~ 1, we have fa 23 ,k r2 ~ 0(A) 
and ki n ,k ri ,k tl ~ (9(A ). However, unlike in the zero- 
flow case, the difference between the high momentum 
transmitted modes fa 2 3 and the high momentum re- 
flected modes k r2 is not negligible for A ^ 1, being 
O (A). Hence, the appearance of such modes requires 
a large momentum change in our system. Therefore, we 
expect the conversion of incident modes fa n into trans- 
mitted modes fa 2 and fa 3 and into reflected modes k T2 
to be disfavoured in comparison with the low momentum 
transmission/reflection channel involving k tl and k ri . In 



other words, we expect that the Wronskian condition ( 38 ) 
will include only the low momentum channel, i.e. 



\Ri\ 



1 



(41) 



From this relation for the reflection coefficient, we 
would conclude that superradiance occurs for all frequen- 
cies in region II (w cr it < to < e<J> ) when A> 1. Note 
that this conclusion certainly does not hold in the case of 
a general dispersive parameter. If the condition A 3> 1 is 
not satisfied, there can be a mixture between the high and 
low momentum channels. Consequently, as discussed be- 



fore, the total reflection coefficient given by equation ( 38 ) 



10 



is not necessarily larger than one. In such a case, a defi- 
nite answer about superradiance can only be obtained by 
solving the differential equations at every spatial point x. 

Having analyzed the case of Wq < 1, let us now fix 
Wq > 1 and assume < u < e&o- As discussed previ- 
ously, there arc two transmission and two reflection chan- 
nels available when Wo > 1 and the scattering problem is 
now described by equation (36). Since we are interested 



in the regime of large A, we calculate the wavenumber of 
the transmitted modes up to next-to-leading order terms, 



1 (u — e$ 



1 

u> ■ 



Wo 
e$ 



2 (1 + Wo) 4 
1 (w - e$ ) 3 



Wo -I 2 (Wo- l) 4 



A- 2 + 0(A~ 4 ), 

A- 2 + 0(A- 4 ). (42) 



By direct substitution of these expressions into equa- 



tion ( 38 ) , one can straightforwadly determine the reflec- 



tion coefficient in powers of A, 



\R^ 



|i? 2 | 2 = l- 



uj — e<i> 



~(\Ti\ 2 -\T 2 \ 2 



(43) 



plus terms of (A~ 2 ) . Note that this second channel 
T2 is present even in the absence of dispersion, being an 
upstream mode which is swept downstream by a super- 
luminal flow. From the equation above, it also becomes 
evident that the relation between the norms |Ti| and |T 2 | 
of the two transmission channels determines the occur- 
rence or not of superradiance. 

Critical case: An interesting situation to be analyzed 
is the critical case oj = u; cr ;t, which corresponds to the 
boundary between regions I and II in FIG. [6] In this sce- 
nario, the background flow Wq, when expanded in powers 
of A, relates to the critical frequency according to the fol- 
lowing expression, 



Wo 



1 



3 ( 

2 l^crit 



e$o)'A-5+o(A-sj. (44) 



Note that the previous analysis leading to equation (41) 



relied on series expansions (see equation (40)) which are 

not valid when Wo — 1 ~ [hr 9 ^ . Therefore, in order 

to analyze the possibility of superradiance in the critical 
case, we cannot use equation (41); instead, we have to 
start from the original Wronskian relation ([Ml). In the 



critical regime, the dispersion relation (see FIG. [7| has 
three distinct solutions. Two of these solutions, denoted 
by k tl and k t2 , have positive group velocities in the lab 
frame and, therefore, are identified as transmitted modes. 
The other solution, denoted by fco, is a degenerate double 
root and, consequently, has a vanishing group velocity in 
the lab frame. In order to obtain the reflection coefficient 
for the scattering problem, we first expand k tl and k t 
power series in A, 



fc tl = a ^°+^(o;- e $o) l A-|+o(A-|), (45) 
fca = -2(w-e$ )*Af +0(1), 



QJA 




k/A 



FIG. 7. The subluminal dispersion relation in the fluid frame 
at the critical frequency. The blue dashed line represents the 
effective frequency Q = ui C rit — e"I>o — kWo at x — > +00. 



and then substitute the obtained expressions into equa- 



tion ( 38 ) . The final result is given by 



\Ri\ 



\R 2 



e$ f,rr, ,2 



(I 



Tir-3|T 2 H, (46) 



plus terms of order O ^A 3 J . Observe again the impor- 
tance of the relative sign of the norm in the two transmis- 
sion channels. Note that the scattering process converts 
incident modes with wavenumber fej n m 10 + 0(A _1 ) into 
transmitted modes of large wavenumber k t2 ~ 0(Aa), 
which can only be balanced by the high momentum re- 
flected modes k r2 w 0(A). Another possibility is the 



conversion of the incident modes ki n 
transmitted modes of wavenumber fc tl 



0(A- r ) into 
(w - e$ )/2, 

which is comparable to the low momentum k ri channel. 
Because of the high momentum change required by the 
second channel (involving k r2 and k t2 ), we expect the 
first channel (involving k ri and k tl ) to be favoured in 
as our scattering experiment. Since w cr i t < e$ j we there- 



fore deduce that the RHS of equation (46) is always 
greater than one in the limit A » 1. In other words, 
low-frequency waves in the critical regime are superra- 
diantly scattered in our toy-model if small subluminal 
corrections are added to the dispersion relation. 
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B. Superluminal dispersion 

We now turn to superluminal scattering processes, 
which can be quite different compared to subluminal ones 
since there does not exist any notion of a horizon or 
mode-independent blocking region for high momentum 
incident modes (the group velocity dw/dk is unbounded 
as a function of k and only the low frequency modes 
which possess quasilinear dispersion experience a block- 
ing region in such flows). We will follow the standard 
treatment [23j of analogue black holes with superluminal 
dispersion and analyze the transmission of an incoming 
wave from — oo through to +00. 

The relevant dispersion relation in the superluminal 
case is depicted in FIG. [8] Given a superluminal flow 
Wq > 1 (solid red curve in FIG. [8|, there exists an in- 
terval of frequencies < oj < w cr i t (region I in FIG. [8| 
for which only two propagating modes are admitted, one 
right-moving and one left-moving in the lab frame. For 
w cr it < oj < e^o (region II in FIG.[8|, however, there are 
four propagating modes, two transmitted right-movers 
and two left-movers. The third possibility is a sublumi- 
nal flow Wo < 1 with < oj < e$o, for which there 
are always two propagating modes (see the black dashed 
curve in FIG. [§. 

Let us discuss first the cases in which only two propa- 
gating modes are available in the asymptotic limit x —¥ 
+00. Hence, the scattering is produced by an incident 
wave from —00 whose frequency oj satisfies < oj < oj cr i t 
(if W > 1) or < oj < e$ (if W < 1). Note that ex- 
actly one of these two propagating modes is a left-moving 
mode. Imposing the boundary condition that no incom- 
ing mode is allowed at +00, we obtain the solution of 
equation (27) corresponding to the scattering problem, 



Re ik r x 



Tie ifc *i x , 



x 

X 



—00, 
00, 



(47) 



plus exponentially decaying channels which do not con- 
tribute directly to the generalized Wronskian current cal- 
culated at x — > ±00. Note that the wavenumbers ki n and 
k r are the same ones that appear in equation ( 23 1 for the 



superluminal W — case and the wavenumber k t repre- 
sents the only available transmission channel. 

On the other hand, if Wo > 1 and the frequency oj of 
the incident wave satisfies w cr it < oj < e<&o (region II in 
FIG. [8]), then there are, in principle, two extra propagat- 
ing channels available (four in total, as discussed above). 
However, because of the boundary condition imposed at 
x — > +00, only one extra transmission channel has to be 
considered (the other extra channel is always left-moving 



at x — >■ +00). The scattering solution is then given by 



f 



e"'" 1 ™'" + tce"* T ", x — > —00, 
T l& ik H x + T 2 e lk ^ x , x -s- 00, 



(48) 



where k in and k r are again given by equation ( 23 ) and k tl 
and k t2 are the wavenumbers of the transmitted modes. 
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FIG. 8. The superluminal dispersion curve in the lab frame 
(A = 1, e$o = 1-5) at x — > +00 for two different flow veloc- 
ities: Wo > 1 (red solid curve) and Wo < 1 (black dashed 
curve). The intervals of frequency indicated refer to the red 
solid curve: region I (green, < oj < w cr it), region II (blue, 
u; cr it < lu < e$o) and region III (light red, non-superradiant 
region) . The intersections of the horizontal lines with the red 
solid and the black dashed curves indicate transmitted modes 
at that frequency. Note that, when x — > —00, the dispersion 
is described by the green dash-dotted curve of FIG. [4] 



Note that we have once again omitted the exponential de- 
caying mode at x — ¥ — 00 since it does not affect directly 
the generalized Wronskian. 

Using equation ( 16 1 to evaluate the functional X in 



both asymptotic regions, we obtain, similarly to the sub- 
luminal case, the following relation between the reflection 
and transmission coefficients, 



|i?| 2 = l- 



A 



VA 2 + 4oj 2 




(oj - e$ - h n Wo) \T n \* 



(49) 
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where the sum is over one or two transmission chan- 
nels, depending on ui and whether Wq > 1 or Wq < 1. 
Here, v gi and v 92 are the group velocities of the trans- 
mitted modes k tl and k t2 , which are always positive 
by definition. Furthermore, it is possible to show, for 
frequencies < u> < e$o, that the effective frequency 
Q = uj — e$ — kt n Wo is always negative for k tl modes 
and always positive for k t2 modes. Therefore, since only 
the n = 1 transmission channel is available for frequen- 
cies lying in region I of FIG. [HJ we conclude that the 
RHS of equation (49) is greater than 1 and, therefore, 



the scattering is always superradiant. 

The situation for Wq < 1 and < u> < e<I>o is simi- 
lar: only the first transmission channel is available and 
superradiance always occurs. However, for frequencies 
located in region II, we cannot so easily conclude super- 
radiance since the extra transmission channel k*„ con- 



tributes an overall negative factor in equation ( 49 1 . To 



obtain a conclusive answer, one would need to know the 
detailed structure of W(x) and V(x) in the intermediate 
regime and solve the equations not only in the asymptotic 
regions but at every point x. 

Large A approximation: In order to better under- 
stand the scattering of an incident wave whose frequency 
is located in region II of FIG. [8j we shall consider small 
deviations from the non-dispersive limit, i.e. A» 1. In 
such a case, we can expand the two transmission chan- 
nels, kt t and k t2 , in powers of A, 



kt, — 



ku = 



1 + W 
u> — e<i>o 



l( W -e$ )3 



2 (1 + Wo) 4 



1 (w-e$ ) 3 A _2 



Wq - 1 2 (W a - l) 4 



A^ + 0(A- 4 ), 



(50) 
(51) 



and substitute the obtained wavenumbers into equa- 
tion (49) in order to determine the reflection coefficient 



for the scattering, 



uj — e$fi 



- (iT^ 2 ~ \T 2 \ 2 ) +<D(l/A 2 ). (52) 



\R\' = l- 



Since we assume uj < e$o m region II, this reflection 
coefficient is larger than 1 whenever \T2\ < |Ti|. As ex- 
plained above, whether this condition is satisfied or not 
in a general model would depend on the detailed struc- 
ture of W(x) and V(x) in the intermediate regime [25] . 
From a practical point of view, in order to maximize the 
potential for superradiance in an experiment with a su- 
perluminally dispersive medium, one should choose the 
asymptotic flow Wo as small as possible while still be- 
ing superluminal (Wo > 1) as this would minimize the 
size of region II and the extra positive-effective-frequency 
transmission channels therein. 

Choosing the flow as such to maximize the Ti channel 
is also consistent with our intuition that scattering favors 
the channel which most closely matches the momentum 
of the reflected mode; in this case the wavenumber fc tl 
is closer to k r than k t2 is. This prediction is confirmed 



in the step function model for which V(x) = e<i>o6(a:) 
and W(x) = Wo0(a;). In such a case, one can impose 
the appropriate boundary conditions at x — discussed 
in the appendix to obtain the following reflection and 
transmission coefficients, 



e$ +o;(Wo-l) . x . 
R ~ u(Wo + l)-e<S> Q +0[A h 

w(W + l)-e$ V ' 



To 



u(W - 1) 
uj(W + 1) - e$ 



O (A- 1 ) 



(53) 
(54) 
(55) 



We can also show that the coefficient correponding to the 



omitted exponential decaying mode in equation ( 48 ) is of 
order O (A" 2 ). 

Comparing T\ and T 2 above and using the fact that 
Wq > 1, one can see that |Ti| > IT2I at zeroth order in 
A, which implies superradiance and confirms our expec- 
tations. Alternatively, one can verify the occurence of 
superradiance by directly analyzing the reflection coeffi- 
cient R above. It is straightforward to see that |i? 2 | > 1 
at lowest order in A. 

Critical case: Another interesting possibility that we 
now consider in detail is the critical regime uj = w C rit- 
This situation corresponds to the boundary between re- 
gions I and II in FIG.[8]and is depicted, in the fluid frame, 
in FIG. [9] The relation between the background flow Wo 
and the critical frequency w cr i t is given by the following 
expression, 



W 



1 + 2 ( Wc «t 



( Wcrit -e$ ) f A"i +o(A"S) . (56) 



Note that, since Wo - 1 - O A 



the A expan- 
sion (|52|) of the generalized Wronskian obtained previ- 



ously is not valid in the present case (check the denom- 
inators in equation (51)). Consequently, we shall need 



different A expansions in order to obtain an appropriate 
expression for the reflection coefficient. 

Like in the critical subluminal case, the dispersion re- 
lation has three distinct roots: the double root ko (with 
vanishing group velocity in the lab frame), a right-moving 
mode (with negative group velocity in the lab frame) and 
a transmitted mode (with positive group velocity in the 
lab frame) whose wavenumber k t is given by 



uj — e<i>o 3 . , .5 . 

h = — ^ — _ 8 3 



O A"3 



(57) 



Applying as a boundary condition the fact that only right 
moving modes are allowed at +00, we obtain, after sub- 
stituting the relevant quantities into equation ( 49 ) , the 



relation betweem the reflection and transmission coeffi- 
cients for the scattering process, 



\R\ 



1 



O 



(A-l) 



(58) 
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airplanes) can be understood in terms of inertial motion 
superradiance [14] . 



IV. APPLICATIONS: AXISYMMETRIC 
SYSTEMS 



Having analyzed superradiance in simple 1+1- 
dimensional toy models with modified dispersion rela- 
tions, we will now discuss how the ideas presented in 
this paper can be generalized to more realistic situations 
based on analogue models of gravity. Our starting point 
is a general 2+1-dimensional, axisymmctric and irrota- 
tional fluid flow with background velocity v given by 



FIG. 9. The superluminal dispersion relation in the fluid 
frame at the critical frequency. The blue dashed line rep- 
resents the effective frequency Q = uj cr it — c$o — kWo at 
x — > +oo. 



Since w cr ;t < e<£>o, we conclude that the RHS of the equa- 
tion above is always greater than one when A > 1. In 
summary, superradiance is expected to occur in the su- 
perluminal critical case for small deviations from the non- 
dispersive regime. 



C. Inertial motion superradiance 

Throughout this paper, inspired by the usual condi- 
tion for rotational superradiance in the black hole case, 
i.e. oj — m£lh < 0, we analyzed only scattering prob- 
lems in which the frequency of the incident mode sat- 
isfies u> — e$o < 0. However, note that since it is the 



effective frequency Q that appears in equation (16), we 



conclude that the amplification of an incident mode can 
also occur for u — e5>o > given that lu — e$o — kWo < 0. 
This kind of superradiance is indeed possible and was an- 
alyzed in Ref. [14] . where it was dubbed inertial motion 
superradiance and compared with rotational superradi- 
ance. In fact, inertial motion superradiance has long been 
known in the literature as the anomalous Dopplcr effect 
and the condition for negative effective-frequency modes 
is referred to as the Ginzburg- Frank condition [24] ■ Sev- 
eral phenomena in physics, like the Vavilov-Cherenkov 
effect and the Mach cones (which appear in supersonic 



v r (r)r + v<t,{r)4> 



A 



B 



-<t>, 



(59) 



where A and B are constants and (r, (f>) are the usual po- 
lar coordinates. Velocity perturbations <5v of the back- 
ground flow can be conveniently described by a scalar 
field ijj, which relates to <5v through equation Sv = V'0. 
We denote the propagation speed of these perturbations 
by c. The idea of analogue gravity is derived from the 
observation that the differential equation satisfied by the 
perturbations "0, i.e. 



(dt 



V)V- 



c 2 V 2 V> = 0, 



(60) 



can be cast into a Klein-Gordon equation in an effec- 
tively curved spacetime geometry. This connection be- 
tween hydrodynamics and gravity is responsible for many 
important results, see Ref. [2] for a detailed review. 

One of the successes of Unruh's [T] original idea of us- 
ing sound waves to study gravitational phenomena is that 
it can be extended to many other physical systems, like 
gravity waves in open channel flows [25] and density per- 
turbations in Bose-Einstein condensates [26] . An impor- 
tant feature of such systems is that equation ( 60 ) is only 



accurate in certain regimes; at sufficiently small distance 
scales (e.g. wavelengths comparable to the fluid depth 
in open channels), the dispersion relation is not linear 
anymore and equation ( 60 ) has to be replaced bjj^] 



(d t + V • v) (d t + y V) ip + <? V> 



+ ^VV, (61) 



where A is a dispersive parameter and the upper (lower) 
sign corresponds to subluminal (superluminal) disper- 
sion. 

Even though superradiance has been subjected to ex- 
tensive studies in the linear regime of analogue models, 
it has never been analysed before in the context of modi- 
fied dispersion relations, as opposed to Hawking radiation 



4 Note that this equation includes only fourth order corrections. 
A full description of the system will possibly include also higher 
order terms. 
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(see e.g. Ref. [5] and references therein). The remarkable 
fact about the toy models introduced in this paper is 
that they can be used to analyze superradiance in re- 
alistic dispersive analogue models of gravity satisfying 
equation (61 ). 



Firstly, the electromagnetic interaction term e$o a P~ 
pearing in our toy models can be mimicked by the rota- 
tional term mfi in axisymmetric analogue models, where 
?7i is the azimuthal number and f2 is the angular velocity. 
It is interesting to note that this duality is manifest in real 
black holes: both electromagnetic [27J and rotational [UI] 
superradiance are possible. Another essential ingredient 
for the occurrence of superradiance in analogue models of 
gravity is the presence of an event horizon, which allows 
no mode to escape from inside the analogue black hole. 
In our toy model, such behaviour is mimicked by an ap- 
propriate boundary condition imposed in the asymptotic 
limit x — > oo. 

It is also important to address the usefulness of the gen- 
eralized Wronskian ( 35 1 in the context of axisymmetric 



systems. More precisely, we are going to show that, if all 
derivatives with respect to x are replaced by derivatives 



with respect to r, then the Wronskian (35 1, when applied 



to solutions of equation (61 ), is independent of r. Indeed, 
by applying the ansatz ifi = (H(r)/y/r) e™^-^ we are 
able to separate equation (61) and are left with a radial 
equation for H, 

H""(r) + a{r)H"(r) + /3(r)H'(r) + ~/(r)H = 0, (62) 

where the coefficients a(r), /3(r) and 7(7") are given by 

1 - 4m 2 



(3 



2r 2 
4to 2 - 1 



±A 2 (c 2 -tv 2 ), 
± A 2 (p(r) - 



2 2 

cr - VI 



(63) 



±A 2 



13m 2 


4 

m 


2r 4 


+ 7T 


! 2 


P(r) 




2r 



■Q(r) 



The functions P(r) and Q(r) appearing in the coefficients 
above are, up to a factor 
and Q defined in Refs. 

as 



t; 2 ) , the same functions P 
. They can be expressed 



P = 



Q = 



+ 1 



c^d_ 
r dr 

ui — 
r dr 



c 
mB 

rv r 

~J2~ 



2iv r 



mB 



2„2 



m c 
^2 



r 
mB 



(64) 



where the upper (lower) sign corresponds to subluminal 
(superluminal) dispersion. Finally, we note that equa- 
tion ( 62 ) is exactly the same as equation ( 28 1 and, more 



remarkably, that the coefficients a, (3, 7 above satisfy 



conditions (30) and (31). Consequently, the generalized 



Wronskian defined in equation ( 35 1 can also be used in 
the context of axisymmetric analogue models of gravity, 
thus completing the connection between our toy models 
and realistic physical systems. 



V. SUMMARY AND FINAL REMARKS 

We have proposed idealized systems to investigate 
multi superradiant scattering processes that are appli- 
cable to sub and superluminal dispersive fields. Perhaps 
the most important theoretical result obtained is related 
to the simplicity of the analytic expression for the particle 
number currents J n of the degrees of freedom / scattering 
channels n (see equation (16)), 



dw 
±oc dk 



(65) 



Notice that the expression above depends only on the 
amplitude, group velocity and effective frequency of the 
particular scattering channel. Moreover, this result is 
universal to all scattering processes discussed in this pa- 
per. Note also that, in principle, we could have normal- 
ized the modes so that the group velocity in ( 65 ) is ab- 
sorbed into the coefficients A n . Doing that would make 



the conservation equations (see e.g. (19)) look much sim- 
pler. However, their dependence on the dispersive pa- 
rameter A would then also be hidden in these new coef- 
ficients. 

Our findings link to standard scattering processes, al- 
lowing a deeper insight into superradiance. Let us con- 
sider scattering of up to four propagating degrees of free- 
dom with a general scattering potential. There are up to 
four propagating degrees of freedom to the left {a, 6, c, d} 
and four to the right {^4, B, C, D} of the scattering poten- 
tial. As a lesson from the analysis carried out before, one 
needs to be careful when assigning the propagation direc- 
tion of each mode since its group velocity is dependent 
on the particular scattering potential and type of disper- 
sion, see FIG. [TUj Furthermore, as long as the scattering 
potential is real, the total current to the left equals the 
total current to the right of the potential, 



Ja + Jb + Jc + Jd — Ja + Jb + Jc + Jl 



(66) 



The scattering potential is amplifying classical and quan- 
tum field excitations if, at the left side of the potential, 
the total outflux (i.e. the reflected current) is larger then 
the total in-flux, J^al > "'total • There are several scat- 
tering coefficients that can be considered: the reflection 
coefficient in each individual scattering channel, 



Pn = 



jJ, 
J, 



ref 



total 



and the total reflection coefficient, 

V 7 rcf 



-Ptotal — 



Tref 
"'total 

J; 



total 



(67) 



(68) 
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FIG. 10. The scattering processes for sub and superlumi- 
nal dispersive fields, with two low (red solid) and two high 
(blue dashed) momentum degrees on each side of the poten- 
tial, are represented by a general scattering process whose 
particle number currents are given by equation (651. 



A sufficient, but not necessary, condition for superra- 
diance, is to demand no influx from the right and to 
strictly require negative effective frequency for the re- 
maining channels to the right — these conditions arise 
naturally at the event horizon of a rotating black hole. 
A surprising result of our analysis is that the presence of 
extra scattering channels can enhance the amplification 
effect. In FIG. [TTJ we illustrate the scattering diagrams 
for some of the examples discussed in the paper. 
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FIG. 11. Schematic representation of some of the scatter- 
ing processes investigated in the paper. The arrows indicate 
the group velocity propagation. The particle number current 
direction is related to the relative sign between this group 
velocity and the effective frequency Q. This information, to- 
gether with equation ( 66 1, can help determine if superradiance 
occurs or not. 



In addition, the present work can be related to recent 
experimental realizations of analogue black holes, some 
of which have even studied the Hawking emission pro- 
cess. In particular, Ref. [5] exhibits the first detection 
of the classical analogue of Hawking radiation using dis- 
persive gravity waves in an open channel flow. One of 



the most important lessons to be learned from that work 
is that even though vorticity and viscosity effects can- 
not be completely removed from the experimental setup, 
they can be made extremely small. In fact, they can 
be reduced to the point that the results predicted by 
the irrotational and inviscid theory match the results ob- 
tained experimentally with considerable accuracy. Based 
on this fact, together with our present results, one might 
ask whether superradiance occurs (and can be observed) 
in the laboratory using dispersive gravity waves. The 
connection between such system and our analysis can be 
seen directly from the full dispersion relation for gravity 
waves in a fluid of constant depth h (30] . 



co 2 {k) 



gk - 
ghk 2 



tanh (kh) 



all 
P 



O (k 6 ) 



(69) 



where g is the gravitational acceleration and a and p 
correspond, respectively, to the surface tension and the 
density of the fluid. If first-order deviations from the 
shallow water limit (kh <C 1) are considered, we recover 
the quartic dispersion relation analyzed in this paper. 

In other words, non-shallow gravity waves impinging 
on a rotating analogue black hole (e.g. a draining 'bath- 
tub' vortex) satisfy equation ( |61[ ) with c 2 = gh and a 
dispersive parameter A~ 2 = \gh /3 — ah/p\. Based on 
the existent analysis of superradiance of linear fields in 
open channel flows 25 , 215] together with our discussion in 



Section IV we expect superradiance to also be manifest 
for non-shallow gravity waves. This relation between dis- 
persion and superradiance in open channels is currently 
being further investigated by the authors and will be the 
subject of a future work. 

Another class of analogue black holes, which was re- 
cently set up in the laboratory [3] and which might be 
used in the future to produce superradiant scattering 
processes, consists of Bose-Einstein condensates [21 126], 
The analogy with gravity arises when one considers the 
Gross-Pitaevskii equation and uses the Madelung repre- 
sentation of the condensate wave function. If the eikonal 
approximation is used and axisymmetry is assumed, it is 
possible to show that perturbations around a background 
condensate obey a superluminal dispersion relation given 
by 



j2 = 4^ 2 ^fc 2 



h 

2m 



(70) 



where m is the mass of a single boson, a is the scatter- 
ing length and no(r) is the background density. Such 
perturbations are described by equation (61 1 with c 2 = 
Airh 2 n a/m 2 and A = 2m/ H. Based on our work, the 
most obvious conclusion we can draw is to expect su- 
perradiance to be manifest also in low-frequency BEC 
scattering experiments. However, since the dispersion 
relation is superluminal, there is no notion of a mode 



16 



independent blocking region (see section III B ) and, con 



sequently, it is not clear if the appropriate boundary con- 
ditions will be sufficient to guarantee superradiance. Ad- 
ditionally, quantized vortices may be present in such sys- 
tems, restricting the angular momentum to integer mul- 
tiples of h. The physics of these quantized vortices with 
respect to superradiance and instabilities is also unclear 
at this point and more investigation is needed to under- 
stand their role in possible BEC scattering processes. 
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Appendix 

Throughout this paper we solve simple models based 
on step functions for the external potential V(x) — 
e$o6(a;) and for the background velocity W(x) = 
WqQ(x), where Q(x) is the Heaviside function. Since 
this function is characterized by a discontinuity at x = 0, 
it is important to analyze what happens to the wavefunc- 
tion, i.e. the solution to equation ^ or equation (28), at 
x = 0. The same problem arises in ID tunneling prob- 
lems in quantum mechanics when the potential barrier in 



the Schrodinger equation is modelled by a step function. 
In such situations, one has to impose the continuity of 
the wavefunction and its first derivative at the disconti- 
nuity point in order to determine the reflection and the 
transmission coefficients. 

Let us first consider equation |9]), which is valid only 
for zero background flows. Following the standard pro- 
cedure used in quantum mechanics (see e.g. Ref. |31|). 
we integrate the differential equation over a small region 
(— e, +e) around the discontinuity x = and then take 
the limit of the obtained expression as e — > 0. Start- 
ing from equation ^ and repeating this procedure three 
times, one can show that f(x) and its derivatives up to 
third order are all continuous at x = 0. 

We also have to deal with the generalization to non- 
zero flows, given by equation (28 1. The situation now is 



more complicate since the differential equation involves 
derivatives of Q(x) (i.e. delta functions). Repeating the 
procedure described above, one can show that the func- 
tion / and its first order derivative are still continuous 
at x — 0. The continuity of the second and third order 
derivatives, because of the delta functions, now depends 
on 6(0), i.e. it depends on the choice of the flow velocity 
and of the external potential exactly at x = 0. Using 
the half-maximum convention of 9(0) = 1/2, we can fur- 
ther prove that the second order derivative f"(x) is also 
continuous at x = and that the third order derivative 
satisfies the following relation, 



/'"(0+^ 



/"'(0-)=A 2 IU 2 /'(0) 

-iA 2 W (cj-e$ Q )f(0). 
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